Progress of Theoretical Physics Supplement 



1 



QCD Critical Points and Their Associated Soft Modes 

T. Kunihiro 1 , Y. Minami 1 and Z. Zhang 2 

1 Department of Physics, Kyoto University, Sakyo-ku, 606-8502, Kyoto, Japan 
2 Mathematical & Physical Science School, North China Electric Power University, 
Zhuxinzhuang, Changping District, Beijin 102206, China 

The mean-field level calculation shows that the QCD matter can have multiple critical 
points incorporating the color superconductivity under charge neutrality constraint due to 
the repulsive vector interaction; this actually implies that the QCD matter is very soft for a 
simultaneous formation of diquark and chiral condensates coupled with the baryonic density. 
Dynamical density fluctuations are analyzed as possible soft modes around the QCD critical 
point using dissipative relativistic fluid dynamics. It is found that the entropy fluctuation 
solely gets enhanced while the sound modes due to mechanical density fluctuations are 
strongly attenuated around the QCD CP, which may suggest a suppression or even total 
disappearance of Mach cone at the CP. 

§1. Introduction 

A unique feature of the phase diagram of quantum chromodynamics (QCD) is 
the existence of a critical point (CP)P>I2) se6) f or a reviewJ3|). In the first part of 
this talk, however, we shall describe possible variants of the QCD phase diagram 
when the color superconductivity(CSC) is taken into account, focusing on a possible 
important role of the vector-vector interaction between quarks.® - ® The mean-field 
level calculation shows that there can exist multiple critical points in the QCD 
phase diagram when the repulsive vector interaction is taken into account,®'® which 
feature can be enhanced when the charge neutrality constraint is imposed® see also 
8) for a similar result in a different context. The possible emergence of multiple 
critical points actually implies that the QCD matter around the phase boundary of 
the chiral-to-CSC transition is soft for a combined excitation of the baryon density, 
scalar condensate and the diquark excitations. 

Around a critical point of a second-order transition, large fluctuations of physical 
quantities are expected. It has been established® ^ that the QCD CP belongs to the 
same universality class(Z2) as the liquid-gas CP, and hydrodynamic modes coupled 
to conserved quantities such as sound modes are softening modes at the CP: We 
show on the basis of the dissipative relativistic fluid dynamics that the dynamical 
density fluctuations or the Brillouin modes are greatly suppressed around the CP, 
while the energy fluctuations or the Rayleigh modes are enhanced when approaching 
the CP.™ 

We shall conclude the report by mentioning some other examples of phase tran- 
sitions that the QCD matter may undergoes and the respective soft modes, which 
could affect the quark spectrum drastically. 
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§2. QCD critical points; alternatives of QCD phase diagram with 
vector interaction and charge neutrality 

When charge neutrality is imposed to the quark matter, a mismatch of the Fermi 
surfaces of the respective flavors arises unless the chiral limit is takenP^ For such a 
system, it is known that the diquark gap has an abnormal temperature dependence 
that the gap increases as the temperature is raised, and has the maximum at a finite 
temperature!^ This is due to the smearing of the Fermi surface by temperature. On 
the other hand, such a system tends to become unstable due to the color magnetic 
instability P) 

It is known that the repulsive vector interaction^ in the Nambu-Jona- 
Lasinio(NJL) modeP^ 1 low-energy effective model of QCEO^ postpones 

the chiral restoration towards larger chemical potentialP'^'^ We note that the 
renormalization-group analysis 13314^5) th e chiral instanton-anti-instanton molecule 
model 2 ^ and the truncated Dyson- Schwinger model of QCD ^ all support the 
existence of the vector-vector four-quark interaction. When the possible transition 
to a CSC phase is considered in the NJL model, it was found that there can appear 
another critical point in the QCD phase diagram.® The relevance of the repulsive 
vector-vector interaction to the chiral and CSC transitions can be intuitively under- 
stood as followsP^ According to thermodynamics, when two phases I and II are in 
an equilibrium state, their temperatures Ti n") pressures Pi t u and the chemical poten- 
tials m t u are the same: T\ = Tn, Pi = Pn, and pL\ = m\. If the two phases are the 
chirally broken and restored phase with quark masses satisfying Mj > M/j, the last 
equality implies that the chirally restored phase has a higher density than the bro- 
ken phase, because /iiji at vanishing temperature are expressed as /i; = J Mf + p 2 F _ , 

(i = I, II), and hence •pp I < PF ri , where pp i is the Fermi momentum of the i-th 
phase. This means that chiral restoration at finite density is necessarily accompa- 
nied by a density jump to a higher density state with a large Fermi surface, which in 
turn favors the formation of Cooper instability leading to CSC. However, since the 
vector coupling in the zero-th component couples to the quark (or baryon) density 
Pb = (<Z7°<7}> it gives rise to a repulsive energy proportional to the density squared, 
i.e. Gvp%/2 with Gy being the vector coupling. This means that the restored phase 
is disfavored energetically due to the vector coupling, and hence the vector coupling 
weakens and delays the phase transition of the chiral restoration at low tempera- 
tures. Thus one also expects that the vector interaction postpones the formation of 
CSC to higher chemical potentials. 

Zhang, Fukushima and Kunihiro^ found that the positive electric chemical po- 
tential jj, e inherent in the charge neutrality constraint plays a similar role with the 
repulsive vector interaction on the chiral phase transition in a four-quark interaction 
model: In a simple two-flavor NJL model, the same two-critical-point structure as 
found in H]) emerges. In addition, positive \i e also represents the magnitude of dif- 
ference between the Fermi spheres of u and d quarks when taking into account the 
local charge neutrality constraint. For an asymmetric homogeneous system with the 
mismatched Fermi spheres, the energy gap of the Cooper paring between these two 
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flavor quarks can increase with temperature,^ as mentioned before. For two-flavor 
neutral CSC phase, this unconventional thermal behavior of the diquark condensate 
can lead to a competition between chiral condensate and diquark condensate, which 
can be enhanced with increasing temperature. For some model parameters region, 
this abnormal competition induced by p e can result in a phase structure with even 
three critical points.^ Thus one would expect that a simultaneous incorporation of 
the repulsive vector interaction and the electric chemical potential under neutrality 
constraint will weaken more significantly the phase transition from the chiral-broken 
to CSC phase. 

Recently, Zhang and Kunihirc® explored the effect of the repulsive vector-vector 
interaction combined with electric-charge neutrality in /3-equilibrium on the chiral 
phase and CSC phase transitions within both two-flavor and two-plus-one-flavor NJL 
models. For the two-flavor case with u and d quarks, a nonlocal NJL modeP2i'U^i i§ 
adopted in[6|): The vector- vector term that is U(2)xU(2) invariant is chosen as 



-V 



-Gv E [(qW-fnqix)? - (q(xh^nq(x)) 2 ] , (2-1 



where q(x) = J dy 4 f(x — y)ip(y). For other notations, we refer to [6]). We remark 
that our choice of the chiral- invariant vector part involves also the axial vector part, 
and the vector (and axial vector) terms ~ {u^u) 2 + (dj^d) 2 have no flavor mixing 
term. We should mention that the other chiral invariant vector interaction (qj^q) 2 
as adopted in® has a flavor mixing and may cause different density dependence of 
the phase diagram from those given in the present work. The nonlocal interaction 
is adopted to deal with the one loop ultraviolet divergence for the calculation of 
the Meissner mass squared at finite temperature.® We take a Lorentzian-type 
form factorpSJ / 2 (p = |p|) = g(p) = —7^, with a = 10, where f(p) is the 

Fourier transformation of the form factor f(x). Oher three model parameters are 
determined by the vacuum physical quantities of the pion mass M n = 135MeV, pion 
decay constant f n = 92.4MeV and the quark condensate — (uu)^' 3 ~ 250MeV. 

We notice that the net quark chemical potential becomes dynamical as given by 
£t a (p) = pL a — 4Gvp a g(p), where p a being the quark density of the flavor a. It is to 
be noted that the induced chemical potential, —AGvp a g(p) and hence the dynamical 
quark chemical potential p a (p) for u and d quarks can be different from each other 
because the electric-charge neutrality makes pd larger than p u and hence pd > p u - 
Thus the mismatch between the effective chemical potentials of u quark and d quark 
becomes 5p = h(p e — ^Gy{pd — Pu)g(p)), which tells us that the difference in the u 
and d quark densities in turn makes smaller the mismatch of their effective chemical 
potentials, thanks to the vector interaction. This effect is found to play an important 
role for the stability against the color magnetic instability. 6 ' 

In the left panel of Fig. [TJ we show an example of the phase diagram with a finite 
vector coupling under the charge neutrality: We use the abbreviations NG, CSC, 
COE, and NOR refer to the hadronic (Nambu-Goldstone) phase with {q%q^) =c a ^ 
and {{qcYails^tabzq))) = A = 0, the CSC phase with A / and a u4 = 0, the 
coexisting phase with o~ u d 7^ and A 7^ 0, and the normal phase with a u d = A = 0, 
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respectively, although they have exact meanings only in the chiral limit. We note 
that there appear four critical points, which are denoted by E, F, G and H. With an 
increase of the vector coupling, the lower two critical points, G and H disappear, while 
the upper two critical points, F and G, remain in the phase diagram. One should 
mention that the critical point H is located on the border between the stable region 
and the unstable region, while other critical points are free from the chromomagnetic 
instability. 

A detailed analysis^ shows that for the parameter set which gives M(p = 0) = 
367.5 MeV and the standard diquark coupling, five different types of chiral critical 
point structures may exist, and the number of the critical points changes as 1 — > 
2 — > 4 — > 2 — >■ with the vector coupling being increased. 



70 

60 

50 

„40 
o 

3, 30 
H 

20 
10 





NOR 




E \ 


NG 


\ F 

CSC 


G v /G s =0.26 




G D /G g =0.75 


1 \ G . . .. 


M(p=0)=367.5 MeV 


/ COE V"'~ 
/ H > Unstable region 





G v /G s = R 
G D /G S = 0.75 
M(p=0)=367.5 MeV 



360 370 380 
\i [MeV] 




400 425 450 
H [MeV] 



Fig. 1. Left panel: The phase diagram for a model-parameter set with Gv/Gs = 0.26 and 
fixed Gd/Gs = 0.75 . Here Gs (Gd) refers to the scalar (diquark) coupling; Gs and other 
parameters are set to gives M(p — 0) = 367.5 MeV. The unstable region with chromomagnetic 
instability is indicated by the dash-dotted curve. There appear four critical points being denoted 
by E, F, G and H. Right panel: The G v dependence of the the unstable region, which tends 
to shrink and move to the high density region as Gv is increased. Both taken from[6j). 



We should notice here that the repulsive vector interaction also suppresses the 
magnitude of the diquark condensate due to the reduced effective quark chemical 
potential. However, the direct effect of the vector interaction on 5[i is more significant 
than that on A, in particular for finite temperature. Thus, as shown in the right 
panel of Fig[TJ we reach a significant findings that an increased vector coupling tends 
to shrink the unsable region toward the hifh-density and lower-temperature region. 
Needless to say, however, the vector interaction may not totally remove the unstable 
region from the phase diagram and hence other mechanism will be still necessary for 
a thorough cure of the magnetic instability. 

The above results obtained for the two-flavor case is not essentially altered even 
when the strange quark is taken into account,^ as shown in FigJ21 where the so-called 
Kobayashi-Maskawa-'tHooft(KMT) termPU - g a i go incorporated 

Ckmt = -K jdet [i> (1 + 75) ip] + det (1 - 75) j . (2-2) 
We should, however, remark that the number of the critical points is sensitive to 
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Fig. 2. A phase diagrams with four critical points for the two-plus-one-fiavor NJL model with 
Gv/Gs = 0.312 and Gd/Gs, under electric-charge-neutrality. Taken from|6}. 



the ratio Gy/Gs and other parameters; other choice of the ratio and other param- 
eters would lead to different phase structures with a varying number of the critical 
pointsP The message of such results we should have is that the QCD matter around 
the phase boundary is very soft for the density fluctuations combined with the for- 
mation of the chiral and diquark condensates along the critical line when the color 
superconductivity is incorporated. 



§3. Density and energy fluctuations around QCD CP 



Now apart from possible variants of the QCD phase diagram, let us examine 
what would be good signatures or observables that reflect the existence of the QCD 
CP. Motivated by the fact that QCD CP belongs to the same universality class 
as the liquid-gas transition, Minami and KunihirJ^ 1 have analyzed the dynamical 
density fluctuations using various dissipative relativistic fluid dynamic equations 
irrespective of the first or second order ones. The dynamical structure factor (or 
spectral function) for the density fluctuation is calculated to be^ 1 

s nn (k,u) 7-1 2Ir/c 2 , i r r B k 2 



({Sn(k,t = o)) 2 ) 7 u 2 + r 2 k 4 ~f v {oj-c s k) 2 + r£k A 

where Ir = K/(n Q c p ) and T B = 3^(7 - 1) + v{\ +\c 2 s T Q (^ - 2r R a P ) with v\ = 
(C + 477/3) /u>o being the relativistic longitudinal kinetic viscosity. The enthalpy 
density in the equilibrium is denoted by wo, while 7 = c p /c v denotes the ratio of the 
specific heats at constant pressure and volume. We refer to Hip for other notations. 
The spectral function has three peaks at frequencies to = and ui = ±c s /c, which 
corresponds to the entropy fluctuations (Rayleight peak), and mechanically induced 
density fluctuation (Brillouin peaks), respectively. Notice that the pre-factor for the 
spectral function for the sound modes is proportional to I/7, which tends to vanish 
when approaching the CP; c p behaves like £ 2_r ? in terms of the correlation length 
£ and the critical exponent n ~ 0.03. Thus one sees that the mechanical density 
fluctuations are attenuated owing to the divergence of the correlation length £ around 
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the QCD CP. On the other hand, the entropy fluctuation in turn is enhanced and 
tends to makes a single peak around the QCD CP in the dynamical structure factor 
of the density fluctuations. 
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Fig. 3. Left panel: The spectral function at t = (T - T c )/T c = 0.5 and k = 0.1 [1/fm]. The 
solid line represents the results using the Landau equation (energy frame) and Israel-Stewart 
equation. The dashed line represents the result using a fluid dynamic equation in the particle 
frameP^ The strength of the Brillouin peaks becomes small due to the singularity of the ratio 
of specific heats. Right panel: The spectral function at t = 0.1 and k = 0.1 [1/fm]. The 
sound modes die out when approaching the CP, irrespective of the relativistic fluid dynamic 
equations used. Taken from llip . 



Such an attenuation of the mechanical density fluctuations may lead to a sup- 
pression or even total disappearance of Mach cone at the QCD CPP If the Mach 
cone formation is confirmed at some incident energy in relativistic heavy-ion colli- 
sions, possible disappearance or strong suppression of a Mach cone along with the 
lowering of the incident energy can be a signal of the existence of the critical point, 
because it may mean that the created matter should have gone through the critical 
region of the CP. 

Eexplicit calculations with equation of motion which admits the existence of the 
critical point is necessary for confirm the fate of Mach cone formation. To make 
a direct connection with experimental observations, we should analyze the density 
fluctuations with the expanding back ground. 



§4. Summary and concluding remarks 



In the first part, after emphasizing that the repulsive vector interaction should 
exist between the quarkswe, we have shown that there is a possibility that the QCD 
phase diagram may have multiple critical points when the color superconductivity 
and the vector interaction are incorporated in the mean-field level. The message 
of this findings is that the QCD matter in the vicinity of the phase boundary at 
low or moderate temperature is very soft for the formation of diquark and chiral 
condensates combined with the baryonic density, as described by aq c q + bqq + cq^q, 
along the critical line. An analysis of such a possibility may involve a dynamical 
Hartre(-Fock)-Bogoliubov theory in the relativistic kinematics, which should be an 
interesting theoretical challenge. 

It was also shown that the vector interaction as given by Eq. (|2-1|) suppresses the 
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chromomagnetic instability related to asymmetric homogeneous 2CSC phase: With 
increasing vector interaction, the unstable region associated with chromomagnetic 
instability shrinks towards lower temperature and higher chemical potential. That 
means that the vector interaction can at least partially resolve the chromomagnetic 
instability problem. 

The dynamical density fluctuations have been analyzed using relativistic fluid 
dynamics, in which the entropy fluctuations are automatically incorporated. The 
sound modes due to density fluctuations are attenuated, and the Rayleigh peak due 
to the entropy fluctuation in turn gets enhanced around the QCD critical point. 
The attenuation of the sound mode may lead to the suppression or even total dis- 
appearance of Mach cone at the CP. Once the Mach cone formation is confirmed in 
experiments of heavy-ion collisions with a high incident energy, possible disappear- 
ance or strong suppression along with the variation of the initial energy can be a 
signal of the existence of the critical point and the created matter went through the 
critical region of the CP. It is clear that there need further explicit calculations for 
a confirmation of this conjecture. 

The present analysis is made in the fluid dynamical regime. In the very vicinity 
of the critical point, we need an analysis beyond the fluid dynamics to take into 
account the non-linear effects. For this purpose, the mode-mode coupling theory® 
and/or dynamical renormalization group technique® should be applied.®® 1 

There are other QCD phase transitions than those discussed in the present re- 
port, and if the phase transition is of a second order or close to that, there should 
exist specific soft modes, which may be easily thermally excited. Such a bosonic 
excitation with a small mass could in turn affect the quark spectral functions: In the 
high-temperature phase above the critical point T c of the color super conductivity, 
the quark spectral function naturally exhibits a pseudo gap® owing to the coupling 
with the pre-formed diquark fluctuations above T c .® The chiral soft modes® may 
also lead to a complicated quark spectral function with three-peak structure above 
but around the critical temperature of the chiral transition.® 
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